Maximizing Fermi surface multiplicity optimizes superconductivity in iron pnictides 
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We study the condition for optimizing superconductivity in the iron pnictides from the lat- 
tice structure point of view. Studying the band structure of the hypothetical lattice structure of 
LaFeAsO, the hole Fermi surface multiplicity is found to be maximized around the Fe-As-Fe bond 
angle regime where the arsenic atoms form a regular tetrahedron. Superconductivity is optimized 
within this three hole Fermi surface regime, while the stoner factor of the antiferromagnetism has 
an overall tendency of increasing upon decreasing the bond angle. Combining also the effect of the 
varying the Fe-As bond length, we provide a guiding principle for obtaining high Tc. 



The discovery of superconductivity in the iron 
pnictides [l| and its Tc up to 55K|^ has given great im- 
pact to the field of condensed matter physics. To attain 
even higher T^, the investigation on the condition for 
optimizing the superconductivity is highly desired. In 
this context, much attention has been paid to the cor- 
relation between Tc and the lattice structure from the 
early stage. In particular, Lee et al. have shown that Tc 
systematically varies with the Fe-Pn-Fe (Pn=pnictogen) 
bond angle, and takes its maximum around 109 degrees, 
at which the pnictogen atoms form a regular tetrahe- 
dron ("Lee's plot")|3[. Theoretically, we have previously 
explained this lattice structure effect within a spin fluc- 
tuation mediated pairing theory on a five orbital model, 
and pointed out that superconductivity is strongly af- 
fected by the Fermi surface around the wave vector (tTjTt) 
in the unfolded Brillouin zone [4]. This Fermi surface has 
been found to be controlled by the pnictogen height /ip„ 
measured from the iron plane. When the pnictogen is at 
high positions, a hole Fermi surface originating from the 
X"^ — orbitalfH'l appears around (7r,7r), and the spin 
fluctuation arising from the interaction between electron 
and hole Fermi surfaces gives rise to a high Tc s±-wave 
paring, where the gap is fully open but changes sign be- 
tween electron and hole Fermi surfaces as was flrst pro- 
posed in ref.Q. When the pnictogen position is low, 
the X"^ — band sinks below the Fermi level, and a 
low Tc nodal s±-wave paringQ or d-wave paring takes 
place[8|. Thus, the appearance of the X^ — Y^ Fermi 
surface around (tt, tt) (which we call 7) is favorable for 
superconductivity, so that the higher pnictogen height 
results in a higher Tc within this scenario. In reffl, the 
saturation of Tc in LnFeAsO in the bond angle regime 
of less than 109 degreesjl, 0] has been attributed to the 
decrease of the lattice constant a, which is found to be 
unfavorable for superconductivity due to the decrease of 
the density of states 4] . 

However, it is becoming clearer that materials hav- 
in g v ery large hpn of > 1.5 A do not have high 
Tc[10|. Recently, for example, Ca4Al206Fe2As2 and 
Ca4Al206Fe2P2, which are variations of the materials 
having perovskite block layers (21311 systems) [Tlj . have 



been found to exhibit superconductivity of 28. 3K and 
17. IK, respectively [l2|. These materials have very large 
hpn (— I.5A), small lattice constant a (~ 3.7lA), and 
consequently a very small bond angle (^ 102 degrees). 
Such a high pnictogen position should give a robust 7 
Fermi surface around (tt, tt) and thus strong spin fluctu- 
ations, and it seems difficult to understand the low Tc 
even if the effect of the small lattice constant is taken 
into account. In this context, a recent band structure 
calculation of Ca4Al2 06Fe2As2 have found an interest- 
ing feature; one of the hole Fermi surfaces around (0,0) 
is missing, resulting in two hole Fermi surfaces, 13i. This 
is found to be due to the very small bond angle. 

Given this background, in the present Letter, we study 
the condition for optimizing superconductivity in the iron 
pnictides, varying hypothetically the lattice structure of 
LaFeAsO. In varying the bond angle a in a wide range 
while fixing the bond length, the number of hole Fermi 
surfaces changes from two to three as a is decreased (as 
discussed as the effect of the pnictogen height in previ- 
ous studies), but as a is decreased even further, one of 
the hole Fermi surfaces around (0, 0) in the unfolded Bril- 
louin zone disappears just as in Ca4Al206Fe2 As2. Conse- 
quently, the number of hole Fermi surfaces is maximized 
around the bond angle where the arsenic atoms form a 
regular tetrahedron. We note that the disappearance of 
the hole Fermi surface was not considered in our previ- 
ous study Q since such a small bond angle regime was not 
considered there. Applying fluctuation exchange (FLEX) 
methodfl^l to the model, we find that the superconduc- 
tivity is optimized in the bond angle regime where there 
are three hole Fermi surfaces. Interestingly, the tendency 
toward magnetism (the strength of the spin fiuctuation 
at zero energy) is not necessarily correlated with the hole 
Fermi surface multiplicity, and therefore the correlation 
between superconductivity and the strength of the spin 
fiuctuations is complicated. We also vary the Fe-As bond 
length while fixing the bond angle, and find that Tc and 
the bond length is positively correlated mainly due to 
the increase of the density of states without affecting the 
Fermi surface. 

We first calculate the band structure of hypo- 
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FIG. 1. The band structure of hypothetical lattice structures 
of LaFeAsO with (a) varying the Fe-As-Fe bond angle while 
fixing the Fe-As bond length and (b) vice versa, (c) The bond 
angle range where three hole Fermi surfaces coexist is shown 
for LaFeAsO and Ca4Al206Fe2As2. 



thetical lattice structures of LaFeAsO using the 
Quantum-ESPRESSO package [l5j|. and construct five 
band models |8:j from maximally localized Wannier 
functions [16|, where we fix the bond length at its original 
length[l[ and vary the bond angle a(Fig[Tray) as has been 
done for Ca4Al206Fe2As2 by Miyake et a/[l3|. The Fermi 
surface is obtained for the band filling(=number of elec- 
trons per site) of n = 6.1. When the bond angle is large, 
two hole Fermi surfaces, ai and a2 are present around 
the wave vector (0,0). As a decreases, the 7 Fermi sur- 



face appears around (7r,7r), and we now have three hole 
Fermi surfaces. This is what has been noticed as an ef- 
fect of increasing the pnictogen height 0, [13, [H| ■ As we 
decrease a even further, the ai Fermi surface disappears, 
and again there are only two hole Fermi surfaces but in 
this case one around (0,0) and another around (7r,7r). 
So the tendency found in Ca4Al206Fe2As2 holds also for 
LaFeAsO. However, the range of the bond angle in which 
there are three hole Fermi surfaces (three hole Fermi sur- 
face zone in Figlljc)) turns out to be much larger for 
LaFeAsO than for Ca4Al206Fe2As2. We find that the 
three hole Fermi surface zone tends to be small for mate- 
rials with large c-axis length. Another point that should 
be mentioned is that the total band width is almost un- 
changed as the bond angle is varied. 

Having understood the effect of the bond angle on the 
band structure, we next fix a = 108 and vary the bond 
length (FiglUb)). The band width is reduced upon in- 
creasing the bond length as expected, but interestingly, 
the Fermi surface is barely affected. In fact, we find that 
multiplying the bands for different bond length by certain 
factors gives nearly the identical band structures, mean- 
ing only the band width, but not the shape, is affected 
by the bond length. 

We now move on to the FLEX calculation. We con- 
sider the standard multiorbital interactions (intraorbital 
and interorbital repulsion, Hund's coupling, and the pair 
hopping interaction), and apply FLEX. In FLEX, bub- 
ble and ladder type diagrams consisting of renormalized 
Green's functions are summed up to obtain the suscepti- 
bilities, which are used to calculate the self energy. The 
renormalized Green's functions are then determined self- 
consistently from the Dyson's equation. The obtained 
Green's function is plugged into the linearized Eliash- 
berg equation, whose eigenvalue A reaches unity at the 
superconducting transition temperature T = T^- Also, 
in order to investigate the correlation between supercon- 
ductivity and magnetism, we obtain the stoner factor as 
of the antiferromagnetism at the wave vector (tt, 0) in the 
unfolded Brillouin zone, which is defined as the largest 
eigenvalue of the matrix C/xo(k = (tt, 0), ia;„ = 0), where 
U is the interaction and xo is the irreducible susceptibil- 
ity matrices, respectively. This value monitors the ten- 
dency towards stripe type antiferromagnetism and the 
strength of the spin fluctuations at zero energy. Since the 
three dimensionality is not strong in LaFeAsO, we take a 
two dimensional model where we neglect the out-of-plane 
hopping integrals, and take 32 x 32 /c-point meshes and 
4096 Matsubara frequencies. 

As for the electron-electron interaction values, we 
adopt the orbital-dependent interactions as obtained 
from first principles calculation in ref [l^ for LaFeAsO, 
but multiply all of them by a constant reducing factor 
/. The reason for introducing this factor is as follows. 
As has been studied in refsl20l- 22 the FLEX for mod- 



els obtained from LDA calculations tend to overestimate 



3 




0.2 
1 



T=0.005eV 
T=0.01eV ■ 



SDW(7r,0) 




95 100 105 110 115 120 

angle (deg.) 




-4 -2 2 4 

bond length (%) 




bond length 



FIG. 2. (a) (Top) the Eliashberg equation eigenvalue for si- 
wave pairing and (bottom) the stoner factor at (7r,0) against 
the bond angle for temperatures T = 0.005 (solid) and O.OleV 
(dashed). The interaction reduction factor is / = 0.6. (b) The 
eigenvalue of the Eliashberg equation and the stoner factor 
against the bond length. T = 0.005eV and / = 0.5. (c) A 
schematic Tc diagram in the bond angle and length space. 



the effect of the self-energy because LDA already par- 
tially takes into account the effect of the self-energy in 
the exchange-correlation functional. When the electron- 
electron interactions as large as those obtained from first 
principles are adopted in the FLEX calculation, this dou- 
ble counting of the self-energy becomes so large that the 
band structure largely differs from its original one. In 



such a case, the spin fluctuations will develop around 
the wave vector (tt, tt) rather than (tt, 0) in the unfolded 
Brillouin zone, which is in disagreement with the exper- 
iments. In the present study, we therefore introduce the 
factor / so as to reduce the electron-electron interactions. 
We vary / within the range where the spin fluctuations 
develop around (tt, 0) and not (tt, tt) . This approach 
would have problems in accuracy if we try to estimate 
the absolute values of Tc (the approach should underes- 
timate Tc), but instead here we calculate the eigenvalue 
of the Eliashberg equation and the stoner factor at a fixed 
temperature, thereby comparing the relative strength to- 
wards superconducting/magnetic instability among var- 
ious lattice structures. 

In FiglJJa), we plot the eigenvalue of the Eliashberg 
equation A for the s±-wave superconductivity and the 
antiferromagnetic stoner factor as at (7r,0) as functions 
of the bond angle while fixing the bond length at its orig- 
inal value. In the large angle regime, the eigenvalue of 
the Eliashberg equation remains relatively small. As the 
angle is decreased, the 7 Fermi surface around (tt, tt) be- 
comes effective and the eigenvalue becomes large. This 
enhancement of superconductivity is correlated with the 
presence/absence of nodes in the superconducting gap on 
the electron Fermi surface (not shown) as has been stud- 
ied previously [3. ItI. [isl. [j^ . As the angle is decreased fur- 
ther, superconductivity tends to be suppressed, although 
the variation of the eigenvalue is not monotonic. Thus, 
the superconductivity is optimized around 110 degrees, 
which is in agreement with the Lee's plotQ. The fact 
that superconductivity is optimized for maximum Fermi 
surface multiplicity can be considered as natural since 
larger number of Fermi surfaces gives rise to larger num- 
ber of pair scattering channels as shown in Fig|3l[a). We 
stress here that due to the quasi-two-dimensionality, the 
density of states is not affected by the size of the Fermi 
surface, so two small Fermi surfaces is more favorable 
than one large Fermi surface regarding the number of 
pair scattering channels. 

On the other hand, the situation is not so simple since 
the stoner factor at (tt, 0), which measures the spin fluc- 
tuations and thus should also affect Tc, is found to be not 
necessarily correlated with the multiplicity of the Fermi 
surface. In the large angle regime of 115 < a < 120, su- 
perconductivity begins to grow as the bond angle is de- 
creased, but the change in the stoner factor is small. This 
can be interpreted as follows. When the — band 
around (tt, tt) sits somewhat away from he Fermi level, 
the interaction between this band and the electron Fermi 
surfaces gives rise to spin fluctuations with finite energy 
that contributes to the pairing interaction but has small 
contribution to spin fluctuations with nearly zero energy. 
This can be related to the NMR experiments in LaFeAsO, 
which do not find any enhancement in the low energy 
spin fluctuations for the optimally doped systems |25l - 
'28|. In the intermediate angle regime (105 < a < 115), 
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FIG. 3. A schematic figure of (a) the difference in the number 
of pair scattering channels between the cases with two small 
or one large hole Fermi surface(s) (b) the difference between 
small and large Fermi surfaces regarding the pair scattering 
processes mediated by spin fluctuations with a momentum 
spread of AQ. 



where the three hole Fermi surfaces coexist, supercon- 
ductivity and the stoner factor is correlated, but in the 
smaller angle regime (a < 105), where ai Fermi surface 
disappears, there is an complicated correlation between 
the two. In this regime, the effect of the 7 Fermi sur- 
face becomes even larger, overcoming the disappearance 
of the ai Fermi surface and resulting in an increase of the 
stoner factor. This works positively for superconductiv- 
ity, and in fact, A is slightly enhanced in the angle regime 
97 < a < 100 for T = 0.005eV. However, in the smaller 
bond angle regime, as continues to increase, but A is de- 
graded. This is mainly because all the Fermi surfaces, 
including the electron ones, become large in this regime. 
Namely, the relative size of the Fermi surface with respect 
to the spread of the spin fluctuations AQ becomes large, 
so that the pair scattering processes occur only between 
restricted portions of the Fermi surfaces, while when the 
size of the Fermi surfaces is comparable to AQ, the pair 
scattering can occur between the entire regime, as shown 
in FiglS^b)^]. Thus, in the small angle regime, although 
the spin fluctuations at (tt, 0) becomes stronger, super- 
conductivity is degraded due to the disappearance of the 
ai and the enlargement of other Fermi surfaces. 

Finally we turn to the bond length dependence. In 
Fig 13b), we plot the eigenvalue of the Eliashberg equa- 
tion as a function of the bond length. The eigenvalue 
monotonically increases with the increase of the bond 
length [3^. Since the change in the Fermi surface is small 
here, the enhancement of the superconductivity is mainly 
due to the increase of the density of states originating 
from the narrowing of the band width. 

In Figl2Ic), we present a schematic diagram in the 
bond angle and bond length space. High is obtained 
when the bond angle is around 110 deg. (the regular 
tetrahedron angle) so that there are (nearly) three hole 
Fermi surfaces and also when the bond length is large so 
that the band width is narrow (sij. 

To summarize, we have studied the condition for op- 



timizing superconductivity in the iron pnictides from 
the lattice structure point of view. As found in 
Ca4Al2 06Fe2As2[l3|, the band structure of the hypothet- 
ical lattice structure of LaFeAsO also exhibits a disap- 
pearance of one of the hole Fermi surface in the small 
bond angle regime, and therefore, the hole Fermi sur- 
face multiplicity is maximized around the bond angle 
regime where the arsenic atoms form a regular tetrahe- 
dron. Superconductivity is optimized within this three 
hole Fermi surface regime, while the (tt, 0) stoner factor 
has an overall increasing tendency upon decreasing the 
bond angle. Consequently, the correlation between su- 
perconductivity and the spin fluctuations becomes rather 
complicated. Combining also the effect of varying the 
Fe-As bond length, we have provided a schematic Tc di- 
agram, which may give a guiding principle for obtaining 
higher Tc materials. 
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